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The phase diagram of the XXZ spin- 1/2 magnet, equivalent to hard-core bosons, under a staggered magnetic 
field and on the Cairo pentagonal lattice is computed at zero and finite temperature by using a cluster mean field 
theory and a stochastic series expansion quantum Monte Carlo. The complex connectivity and the frustration 
lead to unconventional phases such as a 1 /3-ferrimagnetic plateau stabilized by quantum fluctuations as well as 
a 5/12-topological phase induced by a local ice-rule constraint. We also report the presence of a ferrimagnetic 
superfluid and its thermal melting. Finally, we discuss the ferro- and antiferro- magnetic (hopping sign) cases. 



These last years, important discoveries in strongly corre- 
lated physics have been reported in systems where frustra- 
tion plays a central role. Geometrical frustration is very in- 
teresting in that respect. A large variety of unusual phases 
is encountered, from (insulating) exotic spin liquids' to 
superconductivity^. Recent developments on frustrated opti- 
cal lattices of cold atoms open new directions for stabilizing 
exotic bosonic phases such as supersolids^ or Bose-metals^. 

A case of interest is that of ice-rule systems for which 
highly degenerate classical ground states (GS) are governed 
by a local constraint^. For spin-1/2 (boson), it corresponds to a 
fixed number of up (occupied) and down (empty) spins (sites) 
on each elementary brick : a tetrahedron on pyrochlore and 
checkerboard lattice JStll] or a triangle on the kagome^l They 
provide striking features such as charge fractionalization'^. 
Coulomb-gas phases and even magnetic monopoles^'^ . Their 
quantum melting is of broad interest but remains a highly non- 
trivial question; aside from the strength of the interactions, 
the lattice geometry can play a relevant role. Some works 
in this direction revealed exotic phases, e.g. a commensurate 
Resonating Valence Bond supersolid on the checkerboard lat- 
tic^^. Hence, it is natural to search for ice-rule systems with 
more complex geometries such as inequivalent site lattices"'. 

Recently, Ressouche et al. rendered accessible the two 
dimensional Cairo pentagonal structure to experiments by 
proposing the iron-based compound Bi2Fe40g'^. So far it is 
the only known compound with such geometry, a spin-5/2 an- 
tiferromagnet made of identical non-regular pentagons with a 
site dependent connectivity q of 3 (Z3 site) and 4 (Z4 site) 
neighbors as depicted in Fig{T]-a. Note that this Cairo lattice 
is the dual of the Shastry-Sutherland lattice for which the GS 
could be a spin-liquid'^. From the quantum side, substituting 
the iron atoms by copper ones in the Bi2Fe409"^, or creat- 
ing complex optical lattices, the square-based cairo lattice for 
example (Fig[T]-b), is undoubtably a major challenge. How- 
ever, in view of the recent original phases reported in frus- 
trated systems^"^ inequivalent-site structures possess all the 
ingredients to expect unconventional physics. 

In this paper, we study a spin-1/2 magnet under a stag- 
gered magnetic field, or equivalently the extended hard-core 
boson Hubbard model, on the Cairo pentagonal lattice. We 
report a rich phase diagram obtained both at zero and at low 
temperature. We focus on the different insulating phases: 
a topological ice-rule of two bosons per pentagon, a 2/3- 
checkerboard and a pure quantum 1/ 3-ferrimagnetic phase 
with no local constraint. Moreover, a large region of ferrimag- 
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FIG. 1: (Color online), a) The Cairo pentagonal lattice and its 6- 
site unit-cell (shaded plaquette) containing 2 24 sites (dark points) 
and 4 23 sites (light points). The continuously deformed square- 
based Cairo lattice used in this work is displayed in the background, 
b) One of the authorized 5/12 topological ice-rule configurations 
with two spins up (bosons) per pentagon (disks). Thick line: typical 
winding-loop along which a shift of all encountered bosons of one 
lattice site preserves the ice-rule constraint. Shaded plaquette: the 
12-site cluster (2 unit-cells) used in the CMFT. 

netic superfluid is identified as well as its Kosterlitz-Thouless 
(KT) transitior|23 at strong repulsion induced by thermal fluc- 
tuations. We also compare the ferro- and antiferro- magnetic 
cases (sign of the hopping) and discuss the case of the uniform 
magnetic field (nature of the chemical potential). 

I. MODEL AND METHODS 

Spin-1/2 on the pentagonal lattice can be described by an 
extended hard-core boson Hubbard model, with a repulsive 
nearest neighbor interaction V, a hopping t and a chemical 
potential /ij, given by: 

n = -t^{blbj + h.c.) + V^n,nj -^fi,n, (1) 

(ij) (id) i 

where i is the site index, fe| is the creation operator and n,; the 
number of bosons. The coiTespondence is done by the map- 
ping S\ = &| and SI — nii — 1/2 and Eq.|l|) is equivalent 
to a XXZ spin-1/2 magnet with spin couplings = V and 
Jj^ = — 2t, under an effective magnetic field hi. The chemi- 
cal potential /ij is an adjustable parameter^"; this leads to two 
important cases for systems with anisotropic c^, it can be ei- 
ther (a) site dependent or (b) constant. In case (a), if pi is set 
to /i + CiV/2, we have the well-known spin-1/2 magnet under 
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a uniform magnetic fielcH In case (b), the system experi- 
ences a staggered magnetic field hi = fji — CiV/2 depending 
explicitly on the connectivity. Experimentally, this can occur 
in systems with alternating crystal structure and the possible 
presence of a Dzyaloshinskii-Moriya interaction^^. These two 
cases are of great importance in quantum magnetism, but the 
latter case can obviously provide unexpected behaviors as has 
been recently shown in ID materials^''. In the present work, 
we hence focus on constant /i and draw out the rich phase di- 
agram induced by the complex connectivity q . We use equiv- 
alently bosonic or spin language when appropriate. 

To simplify the analysis, we map the parameters on a sphere 
described by two angles 9 and (f) in such a way that p. = sin 0, 
t = cos cos and V = sin0cos(/), where {t,V,Jl) are di- 
mensionless. We also consider a deformed square-based ver- 
sion of the original lattice, see Fig{T]-b. The phase diagram 
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FIG. 2: (Color online). Schematic representation of the insulat- 
ing phases of Fig[3] Dark (light) circles stand for zn, (23) sites 
and the size is proportional to the magnetization (boson density), 
(a) A typical 5/12-topological ice-rule configuration, (b) the 2/3- 
checkerboard and (c) the quantum 1 /3-ferrimagnetic. 

(Fig|3]l is obtained by using two numerical methods: a Cluster 
Mean Field Theory (CMFT'/^ and a Stochastic Series Expan- 
sion (SSE) Quantum Monte-Carlo (QMC)^** at respectively 
zero and finite temperature. The basis of the CMFT is to con- 
sider a finite cluster for which internal bonds of Eq. ([T]) are 
treated exactly whereas the boundary conditions are coupled 
to an external bath. Here, we use a 12-site cluster shown in 
Fig{T]-b (shaded region) and corresponding to two unit-cells. 
The system is diagonalized and solved self-consistently. Note 
that for the triangular lattice, this method gives an excellent 
agreement with the QMC result^. Here, it also allows us to 
consider the frustrated i < case (antiferromagnet) in addi- 
tion to the r=0 properties. The SSE algorithm provides unbi- 
ased quantum simulations for very large system sizes, in our 
case N = l^/ijl x l^jijlwi^h I up to 72 (3888 sites), and 
at finite temperat ure T . Usually, = 21 is enough to focus 
on GS properties^°^'. At very large repulsion however, the 
Kosterlitz-Thouless temperature Tkt significantly drops down 
and a finite temperature transition at small T is expected. In 
this work, the phase diagram is computed at T^^ = 100 > 21 
(up to / = 48) and thermal effects are considered. 

We consider four quantities: x the average number of 
bosons (spin magnetization), ps the superfluid density (spin 
stiffness) implemented via the winding numbers in the SSEPl 
Pn the average number of pentagons with exactly n up-spins 
and the order parameter A/^(q) — {ipo\n{—q)n{q)\'ipo)/N 
with 7i(q) = J2i e^'^ '^'ni performed separately on all sites of 
the non-Bravais square-based lattice a t q = (tt, tt) (Mau) and 
on the Zi sublattice at q = (0,0) (Mx^jSEIj 'pjjg square-based 



Cairo lattice being a depleted square lattice with extra bonds 
(longer bonds on Fig{T]-b), a finite A/(q) is expected even for 
a disordered phase, as entirely explained in ||3T| . 



II. OVERVIEW OF THE PHASE DIAGRAM 

In Fig|3]are depicted the zero (CMFT - dashed lines) and fi- 
nite (SSE - symbols) temperature phase diagrams in the large 
repulsion limit 9/tt > 0.3 where insulating phases appear To 
characterize the different phases, we have considered two cut 
lines respectively at fixed 9 and shown in Fig[3] Since the 
phase diagram is very rich, we briefly introduce it in this para- 
graph before giving more details in the rest of the paper Close 
to 9 = tt/2 [V/t — > 00), only reachable by the mean-field, the 
frustration leads to two magnetization plateaus at a: = 2/3 and 
5/12 (see Figj2]) at T = 0. Surprisingly, when quantum fluc- 
tuations are turned on (9 < tt/2) a third insulating plateau at 
X = 1/3 arises. This insulator is not stabilized at the classical 
limit, is fully driven by the quantum fluctuations and stabi- 
lized by the frustration (t < 0). The quantum melting of these 
lobes leads, in spin language, to a ferrimagnetic superfluid 
(SF) corresponding to different magnetizations (boson den- 
sity) on each sublattice. Uncompensated phases have already 
been reported in systems with complex coordination'^. The 
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FIG. 3: (Color online). Phase diagram in the (6, cj>) plane and the 
comparison between zero (CMFT, dashed lines) and finite (SSE for 
6/-K < 0.485, symbols) temperature. Insulating phases (lobes) are 
enhanced at zero temperature and by frustration (t < 0). The fer- 
rimagnetic superfluid (SF) has a Kosterlitz-Thouless (KT) transition 
to a disordered region (Dis.) at finite T. Lines 1 and 2 are the scans 
used respectively in Fig|4]and Fig|5] 

finite T phase diagram is computed via the SSE method up to 
9/n — 0.485 (circles in Figj3|l. Stronger interaction results 
(shaded regions) are extrapolated. The main difference with 
the T = case is the presence of disordered regions (Dis.) 
due to thermal fluctuations, a Kostertlitz-Thouless transition 
of the SF phases and/or the melting of the insulating lobes. 
Note that small discrepancies between the methods cannot be 
avoided. In the following, we detail the phases of Figj3] In 
particular, we describe the ferrimagnetic character of the SF 
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phase and provide a temperature analysis at strong interaction 
before presenting the unconventional insulators. 



III. SUPERFLUIDITY AND THERMAL FLUCTUATIONS 

As mentioned above, at zero T and for 6 small enough, the 
system is a superfluid {ps 7^ 0) with an on-site magnetization 
(boson density) depending on the connectivity q. We refer 
to this phase as the ferrimagnetic SF^*. We have computed 
the four physical quantities following two cut-lines depicted 
in Fig[3]at both fixed 6 (line 1, Fig|4]i and (j) (line 2, Fig|5]l. On 
Fig|4]are displayed from the top to the bottom x, ps, M{q) 
and p„, as a function of and at finite T. More information 
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FIG. 4: (Color online). Scan along the line 1 of Fig|3]at strong repul- 
sion {Q — 0.4757r), for the cluster size I = 40 and at finite T. From 
top to bottom: the magnetization x (boson density); the spin stiffness 
p3 (superfluid density); the order parameter M'^{q) at (tt, tt) for i\faii 
(all sites) and (0, 0) for M^. (zi sublattice); the average number of 
pentagons p„ with n up spins (bosons). Inset: sublattices used for x 
and (g) indicated by arrows. 



can be obtained by computing x and M{q) on the sublattices 
Z3 and Z4 as well. For (p < 0.27r, Xz^ 7^ Xz^ while ps ^ 



hence corroborating the presence of the ferrimagnetic SF. At 
T = and for 6 large enough (> 0.457r), either transitions 
of first order between two lobes or second order with a super- 
fluid are observed (reinforced at t < 0). At finite T (SSE), 
disordered regions (Dis. in FigjSjl with a finite compressibil- 
ity dx/d<j) 7^ and ps ^ emerge, related to a Kosterlitz- 
Thouless (KT) transition; a finite ps (broken U(l) symmetry) 
in 2D is indeed allowed up to a Tkt temperature. In Fig|5]is 
depicted the size-scaling up to ^ — 72 (upper panel) and the T 
dependence (lower panel) of ps for (f) — 0.1 tt (line 2). From 
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FIG. 5: (Color online). Scan along the line 2 of Fig|3]at ^/tt = O.lof 
Ps as a function of the size / at T = 0.01 (upper) and as a function 
of the T for sizes up to / — 16, 24, 32 and 40 (lower, symbols). 
The dashed line is the 2T/-K line for which a crossing with psiT) 
indicates the Kosterlitz-Thouless transition. Shadings: disordered 
phases of Fig[3] 

the upper panel, a transition between p = and 7^ occurs at 
the thermodynamic limit (TDL). For 2D systems, the KT tran- 
sition is located by the universal jump at P5(7kt) = "^Ty^y/tt 
plus some logarithmic finite size corrections^^. This is cal- 
culated in Fig js] for 3 representative values of 6/tt along line 
2: (i) deep in the SF phase (0.435), (ii) close to the transition 
(0.45) and (iii) in the = region (0.475). We obtain (i) 
T < Tkt, (ii) Tkt 0.0125 close to T and (iii) Tkt < T 
with a KT transition at 9/tt ~ 0.46. The finiteness of Tkt is 
non-trivial and beyond the scope of this paper However, it is 
related to the propagation of defects while doping an insula- 
tor For example, adding a boson in the 5/12 ice-rule phase 
results in creating two defects with zero energy dynamics. 
This would be compatible with a small but finite Tkt at the 
TDL. 



IV. THE INCOMPRESSIBLE PHASES 

The 2/3 -checkerboard phase - On Figj4] at (p/n > 0.3, all 
the pentagons carry 3 bosons (p^ = 1), ps = and M (q) is 
finite for the Z3 sublattice while zero on the Z4 sites. This is in 
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agreement with filled Xz^ sites and empty x^,. This order is 
a simple checkerboard crystal (Fig|2]-b) and is the largest lobe 
of Figj3] with 6'n,ax — O.SStt at </> ~ 0.397r since adding one 
particle costs the energy AV . All the Z3 sites are up spins. 

The 1/3-ferrimagnetic phase - The 1/3 plateau is one of 
the unexpected phases obtained in this geometry which arises 
only under quantum fluctuations (see Figj3]l, both at zero 
and finite T. It is stabilized either under a staggered or a 
constant^'* magnetic field and its expansion is 10 times larger 
(in unit of V/t) for i < 0. It is insulating {ps = 0) with 
no ice-rule constraint (pi 2 7^ 0, Fig|4ji and no broken lat- 
tice symmetry The internal unit cell densities x^.^^ and Xz^ 
mismatch (Fig|4]i showing a ferrimagnetic character at mag- 
netization 1/3 (Fig|2]-c). As displayed in Fig|4]at finite T, the 
1/3 phase is surrounded by two phase transitions, a first or- 
der with the SF for ^/tt ~ 0.01 and a continuous transition 
with the disordered phase. At T = however, the disordered 
phase is not present and a first order transition between the 
1/3 phase to the 5/12 one is obtained (dashed lines in Figj3|l. 
The tip of the 1/3 lobe in Figjsjis located at 6'/7r ~ 0.455 
precisely where the first order transition vanishes. 

The topological ice-rule phase - The 5/12 plateau is sta- 
bilized when the magnetic field is staggered (constant /ij. 
Stable at the Ising limit (CMFT), this phase is robust against 
both quantum and thermal fluctuations, specifically at i < 
(Fig|3]l. On Fig|4j ps = and an average of 2 bosons on Z3 
sites against 0.5 on Z4 ones per unit-cell is observed. M{q) 
is larger than the depletion induced internal structure sig- 
nal, indicating a clear difference between the sublattices^ . 
With p2 — 1, we deduce the presence of an ice-rule of two 
bosons per pentagon. A typical configuration of this 5/12 
phase is given in Fig|2]-a. By labelling the Z4 sites inside a 
unit-cell as and Z4 (inset of FigHi, we identify a finite 
distribution of and x^b w.r.t. (f) related to the degener- 
acy of the GS. Classical zero-energy configurations in stan- 
dard ice-rule systems can generally be connected by quan- 
tum tunneling of a finite number of particles on a closed path. 
This leads, through perturbation theory, to a quantum effec- 
tive Hamiltoni an e. g. quantum dimer models (QDM^fSEEl 
or loop model^'^'-^^'. Here, no such local moves are available, 
only winding-\oo^s invoking the boundary conditions instead 
(see Fig|T]-b). Bosons on such a loop have only two possi- 
ble positions that respect the ice-rule constraint and a tun- 
neling from one to the other results in a new 5/12 config- 
uration. For a cluster of size /, there are Z/4 distinct con- 



tours in each direction {x and y ). For a given configura- 
tion, the number of bosons on such a contour is a conserved 
quantity and each set of these quantities defines a topological 
sector Only the global shift of the bosons along a winding- 
loop can change this number and thus the topological sector; 
the system is protected from local disorder. Starting from 
the most symmetrical 5/12 configuration (Fig{T]-b), the total 
number of winding loops is ^/2. A shift of the bosons along 
one direction cancels the possible winding-loops along the 
other The number of connected configurations is then simply 
O = 2 X 2'-/^. The zero-temperature entropy per site hence 
scales as S /N ~ '''^g;!"^^ and vanishes at the I — )■ 00 limit. 
Since all the configurations are frozen, the phase transition 
is not smooth, as confirmed by the sudden appearance of the 
distribution. We estimate the width of the 5/12 plateau by the 
expansion of this distribution, e.g. 0.13(1) < (p/n < 0.26(1) 
in Fig|4] To our knowledge, isolated sectors have always been 
reported in systems where local moves were also availabld^. 



V. CONCLUDING REMARKS 

We report the phase diagram of spin-1/2 magnets (hard 
core bosons) on the Cairo pentagonal lattice, at zero (CMFT) 
and finite (SSE) temperature. The anisotropic connectivity 
leads, at constant p,i, to a staggered magnetic field^^ Vari- 
ous insulating phases are identified among which a pure quan- 
tum induced 1/3-ferrimagnetic phase, not stabilized in the 
Ising limit {6 — tt/T'/'^. An original 5/12-topological ice- 
rule phase is evidenced, emphasizing the main difference be- 
tween spin-1/2 systems under staggered and uniform fields. 
In this system, the effect of the frustration (t < 0) enhances 
the insulating phases. Finally, a zero vs. finite T comparison 
reveals a KT transition located at strong repulsion as well as a 
partial melting of the lobes. The two methods employed here 
are in good agreement. Open issues remain, such as the com- 
plete description of the phase transitions and the temperature 
properties. Nevertheless, the spin-1/2 Cairo magnet is a very 
promising candidate for exploring new states of matter 
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